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M. ZUPARIC 

Abstract. We show that the scalar product of the phase model on a finite 
rectangular lattice is a (restricted) r-function of the 2-Toda hierarchy. Us- 
ing this equivalence we then show that the wave-functions of the hierarchy 
correspond to certain classes of boundary correlation functions of the model. 



0. Introduction 

In [1], it was observed that the N x N domain wall partition function, Z^, 
of the six vertex model is, up to a multiplicative factor, a t- function of the KP 
hierarchy. In [2], the length M XXZ spin-i chain and its associated scalar product, 
({A}|{/i}), were considered. Restricting either the initial or the final state to a 
Bethe eigenstate, the resulting expression is again a KP r-function. 

In this work we further extend the known correspondences between integrable 
quantum lattice models and classical hierarchies of non linear partial differential 
equations in the following wajQ We show that the scalar product of the phase 
model [4,5], up to a multiplicative factor, is a (restricted) r-function of the 2-Toda 
hierarchy [6,7], where the Toda time variables are power sums of the rapidities. We 
then consider the two types of wave-functions from the Toda theory, and 'w^^\ 
and show that they correspond to a certain class of boundary correlation function 
from the phase model perspective. We additionally give a single determinant form 
for each of these correlation functions. 

In section 1, we recall known results about the finite 2-Toda hierarchy, its con- 
struction from the wave-matrix initial value problem and the r-function as a finite 
bilinear sum of character/ Schur polynomials. In 2, we introduce the phase model 
and the aforementioned lattice model/hierarchy correspondence. Additionally we 
use known combinatorial bijections to express the state vectors of the model as 
weighted sums of various objects. In 3, we use the weighted sum expressions of the 
state vectors to show that the Toda wave-functions correspond to specific classes 
of boundary correlation functions and give their single determinant form. In 4, we 
offer some remarks. 

1. The finite 2-Toda hierarchy 

1.1. Definition of the hierarchy. The details of this section can mostly be found 
in [6,7]. We begin by giving the definition of the 2-Toda hierarchy, with two sets of 
(n — TO — 1) time variables, in terms of four distinct Lax type systems of first order 
differential equations. 

Defining the following shift matrices, 

^lvi,n) = iSk±j,l)k,le{m,....n-l} 
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the general matrix A ^ gl(n — m) is written in the form, 

n — m — 1 
J — — n+m+1 

where m < s < n — 1 denotes the row of the matrix A. A is said to be a strictly 
lower triangular, A = {A)_, if aj{s) = for j > and upper triangular, A = {A)_^_, 
if aj (s) = for j < 0, 

71 — m— 1 —1 

(^)+= E a,(.)A^^,„) , = E a,(.)A^^,„) 
We define two sets of time flows x and y as, 

X = {X1,X2, . . .,Xn—m-l} , V = {^1, J/2, ■ ■ . ,J/n-m-l} 

and introduce L{x,y), M{x,y), Bk{x,y),Ck{x, y) e gl{n — m) where, 
1 

L= J2 b^i''^'y)H^,ny b^{s) = l, Bk = {L'^)+ 

j = — n-\-m-\-\. 
n — m—1 

j=-l 

We deflne the 2-Toda hierarchy as the following Lax type system of differential 
equations, 

(1) d^,L = [Bk,L] , dy,L = [Ck,L] , d^,M = [Bk,M] , dy,M = [Ck,M] 

or equivalently (theorem 1.1 of [6]), the Zakharov-Shabat system, 

dxjBk — dx^Bj + [Bk,Bj] = , dy-Ck — dyj^Cj + [Ck,Cj\ = 

dy.Bk-dx^Cj + [Bk,C.i] =0 

Compatibility conditions. It can be shown that the above systems axe the 
compatibility conditions of the linear operator equations, 

(3) 9,^-1^^°°/"^ = B,-W^(°°/°' , = CjW^°°'°^ 

where W^°°/'^^ = W'^°°/'^\x^y) G GL{n — m) are referred to as wave-matrices. 

1.2. The initial value problem. By defining the constant matrix A G GL(n — 
m) = (aij)i,j=m,...,n-i, such that det [aij]ij=rn...,s-i ^ 0, m < s < n, it is possible 
to flnd wave-matrices, W^°°^ and PF^^) , such that 

(4) W^°^ = W'-'^^A 
where W^°°^ and W^^^ have the specific form 

'n — m — 1 



(2) 



n — m—1 



W 



(oo) 



m<i,j <ri — 1 



k=l 



{w\'^]{i,x,y)^^ 
(wf\{i,x,y)\ 



(5) 



~(oo) 



1 j = '"^^ 







■"'1°'' (^i 7^ const, j — 



j<0 



The remaining non zero entries of the wave-matrices, and are given by 

(proposition 3.1 of [7]), 

det[aij{x,y)]i^^^„,^,Zk,...,s 



wi-\s,x,y)={--^r 



(6) 



det [aij{x,y)]^^ 



-, for 



,j='m,...,s — l 



< k < s — m 
m < s < n — 1 



0</c<n-s-l 
m < s < n — 1 
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where, 



m < i , J < n - 



.1 = exp 



A exp 



and s refers to the row of the entry. The entries of the inverse of the wave-matrices 



^ and (vK(°°)(x,y) 



are similarly given as 



(7) 



wf\s,x,y) = (-1)'=— ^ 
det [aij(a;,i/)] 



for 



i,j—m, ...,s 

det [aij{x,y)]i=m,...,s-i,s+k 

-, for 



< A: < s — m 
m < s < n — 1 



0<k<n~s-l 
m < s < n — 1 



where s now refers to the column of the entry, as opposed to the row. 

The generalized Lax and Zakharov-Shabat systems. From proposition 3.2 
in [7], the following matrices, 

and Bk — {L''} ^ , Ck — {Af*^} , satisfy the linear operator equations (eq. pi, the 
Zakharov-Shabat equations (eq. |2]) and the Lax equations (eq. [T]) which define the 
2-Toda hierarchy. 

1.3. Tau- function of the 2-Toda hierarchy. The r-function, t(s, x,y),is a func- 
tion of the time parameters x and y and an additional parameter, s, which corre- 
sponds to the row number of W^°°\w''^^ or the column number of (W'^°°^) , 

(VF^^-') ^. The derivatives of T{s,x,y) correspond to the entries of the wave- 
matrices and using this fact, we can express the 2-Toda hierarchy in a single integral 
bilinear form. 

Proposition 1. For the function, 

(8) T{s,x,y) = det{aij{x,y)]^^.,^^^_-^ 

the following four relations hold, 



" (oc) / \ 



(9) 



whe 



Cfc {-9s) T-{s, X, y) 
T{s,x,y) 
Cfc (ds) -r(s + l,x,y) 
t{s + l,x,y) 



Cfc [-dyj t{s + l,x,y) 

_ T{s,x,y) 
Cfc [dy) ris, X, y) 
T{s + l,x,y) 



and the generating function for the one row character polynomial, Ck{x), is 
given by, 

oo f n — m — 1 

(10) ^2''Cfc(5) = exp.^ z'xj 

fc=0 I j = l 

Proof. If the above four relations are true then their weighted summations are 
given by, 
(11) 

y^;^fc^.(°o)/„N ^ t(s,^- e(A),y) ^fc^^^(o)^_^^ _ t(s + l,f,y- e(A)) 

^ ^ Tin. T. lT\ 



fe=o T^ys,x,y) 

n — s — 1 , T-* -*t \ \ -'\ s — m 



fc=0 



fc=0 



T{s,x,y) 

r(s,x, y + £(A)) 
r(s + l,x,y) 



J m<i,j<s — 1 
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where, e^A) = (A, . . . ). By using the methods in proposition 3.4 of [7] we 

explicitly obtain, 

r(s, X =F e{\),y) = det [(1 - AA[„,„))±M(f , y)] „<,^^.<^_i 
r(s,f,27=F e(A)) = det L4(f , j7)(l - AA^_„))'''M 

L J T 

which upon expanding as a polynomial in A we obtain the required summations in 
eq. [11] □ 

2-Toda Bilinear relation. The function T{s,x,y) defined in eq. [s] satisfies the 

following bilinear relation, 

(12) 

dX..^_, ("-^-\ ,^ .1 r(s + l,f,y--?(i))r(.'-l,f',y' + r(i)) 



T{s,x,y) T{s',x',y') 



for general s' ,x,x' ,y,y' . The integration § simply refers to the algebraic op- 
eration of obtaining the coefficient of ^. 

Polynomial expressions of the r-function. Rewriting the shift matrix ex- 
ponentials appropriately. 



m<i,j <n— 1 



exp 



C n — m — l ^ n — m — 1 

exp<^ ^ XiAj„,„) I = ^ Cj(^)A[„_„) = 

I ! = 1 J j=0 

n — m— 1 n — m— 1 

- ^ yi(Af„.„)) U E G(-y)(A^n,n))' = (Cw(-y)) 



and using the repeated application of the Cauchy-Binet identity, the r-function 
becomes, 

(13) T{s,x,y)= A{x}{^]X{x}{x)x{t,}{-~y) 

where {A} and {/x} are partitions contained within the box of dimensions (n — 
g^(s-m)^ ^^y^i^x) is the character polynomial given by, 

(14) X{A}(5) = det [CA.+i->(5)]i<,,,<,_„ 
and, 

1.4. Restricting the time variables. In order to make contact with the phase 
model, it is necessary to restrict the time variables in such a way that the r-function 
becomes an element of the symmetric polynomial ring, 

C{[mi, . . . .Us-m]^'-"" , [Vi, . . . ,Vs-m]^'-'"}. 

In the remainder of this work we shall use the convention that, 

• t{x, y) denotes that the time variables are algebraically independent. 

• t{u, v) denotes that the time variables are algebraically dependent, and 
t(u, v) is an element of the aforementioned symmetric polynomial ring. We 
shall refer to t{u, v) as a restricted r-function. 

Schur polynomials. Changing from time parameters to symmetric power sum^ 

(15) Xk = ipfc(tti, . . . ,Ms-m) , — yk ^ ipfc(i>i, . . . ,Vs-m) , fc e {1, . . . ,n — m — 1} 



^We define the symmetric power sum of order k as, Pfc(«) = 
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the one row character polynomials, Cii^) ^^J^d Ci(~2/*)i become complete homoge- 
neous symmetric polynomials]^ 

Hence the character polynomials in the r-function become Schur polynomials, 

(16) r{s,u,v)^ yl{A}{M}'S'{A}(^)S{^}(iT) 

2. The phase model 

Most of the material in this section can be found in [4, 5] . Consider the bosonic 
algebrsj^ generated by the three operators (j), (f>^ and A'^ that satisfy the following 
commutation relations, 

[4>, 4>^]=n , [N, 0] = -0 , [iV, 0t] = 0t 

where tt — |0)(0| is the vacuum projector. The one dimensional Fock space, F, for 
this algebra is formed from the state \n), where the label n G Z+ U {0} is called an 
occupation number. The action of the operators on elements of the Fock space are 
given by, 

4>'*\n) = In + 1) , (l>\n) = |7i - 1) , N\n) = n\n} 

The action of the operator on the vacuum state, |0), annihilates it. 

We now extend the above bosonic algebra and consider the tensor product, 

F = Fo®Fi®---®Fm 

We introduce the operators, (pj, and Nj, < j < M, that act on ¥j as, 

= /() ® /i ® • • • i8) i8) Ij+i (8 ■ ■ ■ <8 hi 

and similarly for 0j and Nj, where Ik is the identity operator in F^. The commu- 
tation relations are given by 

(17) = TTjSjk , lNj,(t)k] = -(pjSjk , lNj,(t)l] = (f>'l5jk 

where each operator of index j acts on the corresponding indexed Fock vectors, 

{(pj)""''"' |mj)j = \nj)j for < Uj < 

(j>]J \mj)j = \nj).j for Uj > rUj > 

Nj\mj)j = rnj\mj)j 

and (pj annihilates the vacuum state \Q)j. The state vectors, \np)j, and the corre- 
sponding conjugate vectors, (n^lfe, are orthonormal, 

(19) {nr\np)k,j = SprSjk 

2.1. Algebraic Bethe ansatz. We define the phase model through the following 
local i-operator matrix. 



(20) L,{i 



hj (u) 



- 0^ 

u ^3 



Cj(u) dj(u) J \ u 

where w G C. Naturally associated with Lj{u) is the 4x4 i?- matrix, 

/ f{u,v) \ 

g{u,v) 1 

9{u,v) 

V f{u,v) J 



(21) R{u,v) 



■^Thc complete homogeneous symmetric polynomials, hk{u), are generated by X]fcLo ^''^ki''^) ■ 
nj=]" i—\u - ~ I ^ 1 additional information see section 1.2 of [8]. 

"^We note that this algebra is the q = limit of the g-boson algebra. 
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where f{u,v) = ^_ 2 and g{u,v) = -^^ 2 ■ L and R satisfy the following inter- 
twining relation, 

(22) R{u,v)[Lj{u)® Lj{v)] = [Lj{v) ® Lj{u)]R{u,v) 



22 1 the monodromy matrix and 



The monodromy matrix. The monodromy matrix, T{u) , for the phase model is 
introduced as the ordered product of all (M + 1) L-matrices, 

(23) T{u) = Lm{u)Lm-i(u) . . . L,{u) = ^[^j ^[^j 

Using induction on the intertwining relation (eq 
the i?-matrix satisfy an equivalent intertwining relationship, 

(24) R{u,v)[T{u)®T[v)] = [T(v) ® T{u)\R{u,v) 
which generate sixteen non trivial algebraic relationships. 

Non local creation and annihilation operators. Focusing on the operators 
B{u) and C(u), we consider the operator N = J^jLo^j' which measures the to- 
tal occupation number of the state. Applying this operator to B{u) and C{u) we 
obtain, 

(25) NB{u) = B{u) + 1} , NC{u) = C{u) |iV - l| 

Thus the operator B(u) is a creation operator of the phase model, where one appli- 
cation on a state vector increases the total occupation number by one, while C(u) 
is the opposing annihilation operator of the phase model, where one application to 
a state vector decreases the total occupation number by one. We note that C{u) 
annihilates the total vacuum operator. 

Equivalently, the roles of B{u) and C{u) are reversed when applied to the con- 
jugated vacuum vectors. 



2.2. N-particle state vector and its conjugate. We construct the A^-particle 
vector, I^'a/), by repeated application of the construction operator B on the vacuum 
vector, 

\^m) ^ Biui)...B{uNm 

where the total occupation number of |^'m) is N. An alternative form for the 
A'^-particle vector is given by, 

M M 

(26) |*m)= J2 /{no,....nM}Wn('^l)"'0l°>^-= E fw(^)W 

0<iio,iii,...,iiM<" fc=0 j=0 {A}C{M)" 

where the partition {A} is constructed from the occupation number sequence in the 
following manner, 

(27) {A} = (M"*^(M- l)""-i,...,l"\0"'') 

Similarly, we construct the conjugate A^-particle vector, {^m\, by repeated appli- 
cation of the annihilation operator C on the conjugate vacuum vector, 

(*Af| = {0\C{vn)...C{vj) 
where the total occupation number N, and an alternative form is given by, 

M A/ 

(28) {*m|= J2 5{no,...,n„}(^?)0(O|, n('^^r'" = E 9wiv){\\ 

0<no."l. i = fe=0 {A}C(M)« 
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Schur polynomial expansion of state vectors. In [9], Tsilevich derived the 
following Schur polynomial forms for and g{x}iv) in eqs. 26 and 

(29) 



28 



/ N \ / N \ ' 



S{\}{vi ^, ■ ■ • ,"iv^ 



In order to proceed, we need to express the state vectors as weighted sums of lattice 
configurations, plane partitions and semi-standard tableaux. 

2.3. Weighted sums of lattice paths. Consider N non-crossing column strict 
lattice paths on the {M + 1) x 2N lattice, where the paths begin at the first N 
bottom-most horizontal edges, (— iV, 0), (— iV + 1, 0), . . . , (— 1, 0), and end at the 
final N top-most horizontal edges, (1, M), (2, M), . . . , {N, M), respectively. Each 
path is restricted to move either up or right, and no two paths can cross or occupy 
the same vertical edge. We give a typical example in fig. [T] 



(1.4) 



Figure 1. Generic lattice path configuration for Af = 4, A'^ = 3. 

Given an allowable lattice path configuration, we assign each of the four possible 
vertices a letter, as indicated in fig. [2] 



/\ /\ \y \/ 

\/ /\ \/ 
abed 



Figure 2 . The four vertices used to construct lattice paths 



Following section V of [5], it is possible to express the coefficient of the conjugate 
state vector, g{na,...,nM}(^)' following weighted sum of lattice paths. 



(30) 



E 



"n 



allowable paths 
(M + 1) X N latti< 



where the sum is taken over all allowable paths in the (M -I- 1) x A'^ lattice under 
the conditions, 

• Uj^ paths start at {—N, 0), {~N — 1,0),..., {—N + nj^ — 1, 0) and end at 
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• Uj^ paths start at {—N + iij^ , 0), . . . , {—N + "-j! ~ ^^ ^) ^^^^ ^'^^ 

(-1,J2). 

• This procedure continues until we finally have rij^, paths starting at (— + 

Y!iZlnji,Q),...,{-\,Q) and ending at (-l,jfc). 

The powers tf and tf , \ < I < N , are equal to the number of d and a vertices 
respectively in the Ith. column. 

Similarly, the coefficient of the state vector, /{no,...,„„}(u), can be expressed as 
the following weighted sum of lattice paths. 



(31) 



/{n,^,...,n,^}(M) = 



E 



.,N 



allowable paths 
+ 1) X N latti< 



where the sum is taken over all allowable paths in the {M + 1) x lattice under 
the conditions, 

• riji paths start at (l,ji) and end at (1, M), (2, M), . . . , (nj^ , M). 

• rij^ paths start at (1, j2) and end at [rij-^ + 1, M), . . . , (X]/=i '^ii ; ^'^)- 

• This procedure continues until we finally have rij^ paths starting at (1, jfc) 
and ending at {Y!IZi n^, + 1, M), . . . , {N, M). 

2.4. Weighted sums of plane partitions. A plane partition, f^j.k, is an array of 
non negative integers such that, 

T^j,k > T^j + l.k and TTj^k > 7rj,fe + i 

If we restrict the size of the array to be x A^, and restrict the maximum of 
any integer within the plane partition, TTi.j < M, the plane partition is said to be 
contained within a box of side lengths N x N x M . 

A typical example of a plane partition within a box of 3 x 3 x 4 is given by the 
following 

(32) TT*'' > = 3 1 1 

V 2 1 1 / 

Note that the diagonal entries of any plane partition always give a partition in 
usual sensc|^ The graphical representation of a plane partition in a x iV x M 
box is given by considering rhombus tilings of a (A^, A^, M) semiregular hexagon. 
The plane partition, tt^^ ^, is represented by fig. [sj where each representation can 
be generally constructed entirely from the three types of rhombi given in fig. [4] 




Figure 3. Graphical representation of the plane partition vr^^ K 



'We shall use the following plane partition, tt^^ J, as a running example in this section. 
'In the above example we obviously have {A'} = (3, 1, 1). 
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Figure 4. The three types of rhombi used to construct plane partitions. 



Correspondence between plane partitions in a N x N x M box and N 
non-crossing column strict lattice paths on the (M + 1) x 2N lattice. The 

jth path of the lattice configuration can be thought of as the jth column of the 
array tt^. As an example, consider the array, tt^"^ \ which is in correspondence with 
the lattice path configuration shown in fig. [T] 

The bottom left entry, Trl\^ — 2, corresponds to the first horizontal section 
of the first path on the second row, and the remaining entries of the column, 
(""2^^ = 3,7r|;'^'^ = 3), correspond to the remaining horizontal sections of the first 
path, on the third row. 

There exists a similar correspondence between the second and third lattice paths, 
and the second and third columns of the array tt^'^ ^ respectively. It should be clear 
how this process is generalized for any N and M. 

Due to the above correspondence, when considering lower diagonal plane par- 
titions in an N X N X M box, (equivalently, the left hand side of the {N, N, M) 
rhombus tiling), we obtain 

(33) E vr'^-vt'K..v'y- 

lower diagonal 
plane partitions 

where the sum is taken over all allowable lower diagonal NxNxM plane partitions, 
(left hand side (A'', A'', M) rhombus tilings) , and the diagonal terms are given by 
the partition representation of the corresponding occupation number sequence (eq. 



27 1. The powers if and If, 1 < I < N , are equal to the number of d and a rhombi 
respectively in the Ith column of the left half rhombus tiling. 

Similarly, when considering upper diagonal plane partitions in an N x N x M 
box, (equivalently, the right hand side of the {N,N,M) rhombus tiling), we obtain 

(34) /{n,,,...,n,,}(^)= E ur'^'u't'K..uf-'« 

upper diagonal 
plane partitions 

where the sum is taken over all allowable upper diagonal N x N x M plane par- 
titions. Again, the diagonal terms are given by the partition representation of the 
corresponding occupation number sequence. 

2.5. Weighted sums of semi-standard tableaux. We now give the final alter- 
native forms of the state vectors using the following correspondences. 

Correspondence between upper diagonal plane partitions, '!t^\ and semi- 
standard tableaux of descending orderj^ T^"^^. We begin by considering a gen- 
eral upper half plane partition, tt^^^, and construct a partition using the diagonal 
entries, 

{A} — {tTi,!, 7r2,2, ■ . ■ , 7i"]V,Jv} 

Considering the next upper diagonal entries, tt^j+i, we construct the skew diagram, 

{Mi} = {^1,1 ^ 1"1,2, 7r2,2 — 7r2,3, . . . , 7r]V-l,JV_i — TTjV-l.JV, 7rjv,jv} 

^Semi-standard tableaux are commonly of ascending numerical order, however, descending 
numerical order is the most convenient convention for the purposes of the next section. 
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and place the integer 1 in the vahd regions of the skew diagrarrj^ We then consider 
the next upper diagonal entries of the array, 7rjj+2, and construct the skew diagram, 
{^2}, 

{^2} = {""l,! — 7ri,3, 712,2 — 712,4, • • • , 7rjv_2,JV-2 — 7rjv,JV-2 , 7rjV-l,JV-l , 7rjv,jv} 

and place the integer 2 in the valid regions of the skew diagram that have not 
already been occupied by previous steps in this process. This process continues 
until the partition contains the numbers {1,...,A^ — 1}. We then fill the remaining 
boxes in the partition with the integer N, thereby constructing a valid descending 

semi-standard tableau T^^^ from the upper diagonal plane partition tt\^^ . 

As an example, consider the array, tt^ , given in the past examples where {A'} = 
(3, 1, 1). The construction of the corresponding descending semi-standard tableau 
is given in fig. [sj 

In (a) we construct the partition {A'} = (3, 1, 1). In (b) we construct the skew 
partition {^1} = (3, 1, 1) — (1, 1, 0) and place the integer 1 in the valid regions of 
{/ii}. The partition (1,1,0) was obtained from the first upper diagonal entries of 
TT^ . In (c) we construct the skew partition {^2} — (3,1,1) — (1,0,0) and place 
the integer 2 in the valid regions of {(12} that contain no integers. The partition 
(1, 0, 0) was obtained from the second upper diagonal entries of . In (d) we place 
the integer 3 in any remaining entries of {A'} that don't already contain integers, 

forming the valid descending semi-standard tableau T^^ ^ from the upper diagonal 
plane partition tt^ ^ . 









1 1 




1 1 


3 


1 1 




— > 

(a) 


1 

1 


— B> 

(b) 


2 
1 


^ 

(c) 


2 

1 


(d) 



Figure 5. Construction of the tableau T 



Thus based on the above correspondence, another valid combinatorial definition for 
/{A} (w) is the following, 

(35) fwW^}^^.! ■■■u^" 

where the summation is over all semi- standard Young tableaux of shape {A}. The 

powers, tj, give the weights of T^^\ which count the number of times j appears in 
the tableau. Note that this expression is in accordance with eq. [29] 

Correspondence between lower diagonal plane partitions, Tr^\ and semi- 
standard tableaux of ascending order, T^'*'^. Using an equivalent algorithm 
as described above, except this time applying a numerically ascending convention, 
we obtain the required correspondence. Using tt^^ ^ as an example yet again, the 
construction of the corresponding ascending semi-standard tableau is given in fig. 

El 

In (a) we construct the partition {A'} — (3, 1, 1). In (b) we construct the skew 
partition {I'l} — (3, 1, 1) — (3, 1,0) and place the integer 3 in the valid regions of 
{I'l}. The partition (3, 1,0) was obtained from the first lower diagonal entries of 
TT^ . In (c) we construct the skew partition {1/2} ~ (3,1,1) — (2,0,0) and place 



'In ascending tableaux, A'^ would be placed instead of 1. 
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the integer 2 in the vahd regions of {^2} that contain no integers. The partition 
(2, 0, 0) was obtained from the second lower diagonal entries of tt'*' . In (d) we place 
the integer 1 in any remaining entries of {A'} that don't already contain integers, 

forming the valid ascending semi-standard tableau t|'^ ^ from the lower diagonal 

plane partition ^ . 




Thus we immediately obtain, 

/Qns ~2ti+M — 2t2 + M -2tj^+M 

(36) fl{A}(w)=2^"i V2 ■■■vj^ 

+ 

where the summation is over all semi-standard Young tableaux of shape {A} of 
ascending numerical order. 

2.6. The scalar product. We now consider the scalar product, S(7V, M|m, 1/), of 
the phase model which is defined as the inner product of the state vectors, 

(37) §(iV,M|u,«) = {*m1*m) 

Using the algebraic expressions from eq. [24] it is possible to obtain the following 

determinant expression, 

(38) 



<m,l<N 



/ N \ 

S(iV,M)= n (^7^) (^;^)det[W-i(-i,«?)], 

Alternatively, considering the Schur polynomial expansion of the state vectors, 

(39) HN,M)=il[^] J2 Si^y{ui,...,u%)Six){v^\...,v],^) 

\j = l ^ I {A}C(M)« 

2.7. Restricting the 2-Toda tau-function to obtain the scalar product. 
Proposition 2. The scalar product of the phase model for general N and M is, up 
to an overall factor of (0^=1 Tt) ' restricted r-function of the 2-Toda hierarchy 
with = 5{x}{^} and s = n—M = m+N , where m and n are free parameters. 

Proof. Beginning with the unrestricted r-function, 

T{s = n-M = m + N,x,y)^ ^ X{a} (k)X{a} (-j7) 

{A}C(M)" 

and performing the following change of variables, 

(40) Xk ^ ipfc (wi, • • • , u%) , -yk ^ ^Pk (i^r^i • • • , "iv^) ,l<k<N + M- l 

we obtain the required result. □ 

The above result only considers one value of s. Let us now consider the family 
of corresponding restricted r-functions for other values of s. We begin by clarifying 
some known facts about the family of unrestricted r-functions. 
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• The entire family is given by {rs=„+i(f , y) , Ts=„+2 (x, y), . . . , rs=„(f, y)}. 

• Different values of s do not change the amount of, (two sets of n — to — 1), 
time variables. 

We now compare this to the case of the family of restricted r-functions. 

• The initial r-function, t(^s — n — M = m + N, {u^} , has two sets 
of -/V + M — 1 time variables, but each set is constructed from N symmetric 
variables. 

• The introduction of the condition s — n — M = m + N means that as s 
changes, so to do M and TV. By considering the change in the dimensions 
of the partition, we can obtain how M and N change with s. 



(41) s ^ s ± r 



N ^ N±l 



• Consequently, although the number of time variables do not change with 
each s value, different values of s do change the amount of symmetric 
variables that the time variables are constructed from. 

An example. Consider the complete family of unrestricted r-functions for n = 5 
and m = 1. In this case each r-function contains two sets of 3 time variables, 
{x,y} {xi,x2,X3,yi,y2,y3}, 

!T'3=2 = ^ 0{A} , rs=3 = ^ 6{A} , rs=4 = ^ 0{A} , rs=5 = 6{0} > 
{A}C{3} {A}C{2,2} {A}C{1,1,1} J 

where 6{a} — X{\}{^)X{\}{~y)- The main question now is, if one r-function in 
a family has been restricted to form a scalar product with a certain M and N 
value, can the remaining r-functions of the family also be restricted to form scalar 
products with valid M and N values? 

Naively performing the corresponding restrictions to the above family of r- 
functions we obtain the following family of scalar products, 

3„/iV = l\ 2 /7V = 2\ 1 /7V = 3\ o^fN = 4: 



where — (11^=1 IT ) ■ This example illustrates an important issue. We re- 
member that each r-function contained within a family must contain the same 
amount of time variables. Furthermore, it is a requirement that these time vari- 
ables be the same for each value of s to ensure that the r-functions obey the bilinear 
relatiorj^ If this is to be the case for the above example, we have the following set 
of equations that must be satisfied, 

(^1) j = Pfe \ ) A^i) ) = ■■■ = Vk \ \u^ ),..., (M4 ) j 
[vx] ^=Vk\(vx) ,(«2) j = ■•■ =Pfe (("1 ) ,...,(^4) j 

for 1 < fc < 3. A simple check will reveal that only the trivial solution exists, 
meaning that all but one of the symmetric variables are set to zero. Thus, at a 
first glance, the answer to the question is no, due to the fact that the t -functions 
in the family all need to contain the same time variables, be they independent or 
restricted. 

We now generalize the above example. 



The r-functions obviously must obey tlic bilinear relation. 
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Proposition 3. The system of equations, < I < M — 1, 

It? H h u%+i = ^i? H \- ^^% 

lit H h u%+i = ^1 H \- fJ.% 



2(]V+M-1) 2(]V+M-1) _ 2(JV+M-1) 2(JV+M-1) 

permits only the trivial solution, i.e. uj,. = for j G {1, . . . , N}, and the remain- 
ing I of the m| 's are equal to zero. 

Proof. We begin by considering the first N + I equations in the system, the 
remaining equations will follow easily. We note that the left hand side of these 
polynomial equations exist in the symmetric polynomial ring C[u?, . . . , 
Consider now the polynomial ring C[si, . . . , sjv+j], and recall that the fundamental 
theorem of symmetric polynomials states that there exists an isomorphism between 
the two rings, C[m?, . . . , UAr+;]'^"+' = C[si, . . . , sat^/], with the isomorphism send- 
inj^ Pj (uf, . . . , uj^_^_i) Sj, j = {1, . . . ,N + I}. Hence the system of N + I 
equations in the isomorphic polynomial ring, C[si, . . . , sat+j], is linear and has one 
solution. 

Thus in the ring C[ui, . . . , , the system contains one base solution, and 

every possible permutation of that base solution (since the polynomial ring is sym- 
metric), leading to a total of {N + 1)1 possible solutions. Since we already trivially 
know {N + ly. solutions to the system, u^^ = /i^ for G {1, . • . ,N}, and u^^ = 
for k {N + 1, . . . , N + I}, this means only the trivial solution exists for the first 
N + I equations. It remains to note that the remaining AI — I — 1 equations are 
solved by the {N + l)\ solutions. □ 

Using the above result the following lemma comes almost automatically. 

Lemma 1. Assume we have a particular family of unrestricted r-functions with 
particular m and n values, 

(42) {r™+i(f,j7),r„+2(f,j7), . . . ,T„{x,y)} 

The process of restricting the entire family so that each T-function corresponds to 
a valid scalar product expression. 



has potentially two (ill) effects. 

• If each of the above scalar product expressions has two sets of N (N is not 
constant for each scalar product) symmetric variables, then the 2 sets of 
n — m — I = N + M — I time variables of the restricted T-functions are no 
longer equal, and therefore the bilinear identity is no longer valid. 

• // we enforce that the time variables be equal, then we only have two sets of 
one symmetric variable for each of the scalar product expressions. 

Arguably both scenarios are pointless, so it makes sense to use the results of propo- 
sition^and only consider restricting one T-function in any family. 

Proof. Applying proposition [2] on all the unrestricted r-functions in eq. |42j 
we instantly arrive to the expression in eq. |43[ Analyzing any two of the above 
scalar product expressions, (with two sets of N symmetric variables), the results of 
proposition [3] state that the symmetric power sums, and hence the time variables, 
cannot be equal. Thus the first point in this lemma becomes obvious. Furthermore, 



"'^''We could use any basis symmetric polynomial, 
of [8] for further details. 
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from proposition [3] the only way for the time variables to be equal is if we trivialize 
the power sums as indicated in point 2 of this lemma. □ 

3. The Toda wave-vectors 

In this section we shall show that restricting the wave- functions of 2- Toda hierar- 
chy give an alternative method to calculating certain classes of correlation functions, 
and thus have a natural combinatorial meaning. In order to proceed we shall first 
give necessary definitions of skew Schur polynomial^"] 

Skew Schur polynomials. Given a set of variables u — . . . ,UAr) and the 
partitions {A}, {/i}, such that {A} I) the skew Schur polynomial, S^y^ f^^:^{u), 
is defined as, 

= det[h\^-fj,-+j-i{ui, . . . , UN)]i<i,j<N 

where the sum is given over all possible (ascending or descending) semi-standard 
skew tableaux of shape {A — /i}, and the tj give the weights of the tableau. 

3.1. Wave- functions - I. Considering the li'"^ class of wave- functions, from eq. 
|9]we obtain, 

r(s)u>^°'(s) = J2 X{A}(5)X{A}/{fc}(-y) 

{A}C(,i-{a + l))((s + l)~™) 

where we have used the following result, 

(45) Cj(-%)x{A}(-y) = X{\}/{j}{-y) 

for all partitions {A} such that {j} C {A}. 

Thus the upper triangular wave-matrix, W^*-*^-*, has entries of the form, 

^'"^-\^) E X{A}(f)X{A}/{.-.}(-37)) 

V ^•^^A}C(„-(, + l))(0 + l)-") /„<,,fe<„-l 

Constructing skew A^-particle conjugate state vectors. Consider the follow- 
ing conjugate state vector, 

mkC{v2)...C{vN)^{H'i^\ 

The allowable partitions of this conjugate vector are given by the following result. 
Proposition 4. 

<*i?^i= E <Af(Al 

{A}C{(M)(™-1) ,1-} 

Proof. Consider the non crossing column strict lattice path interpretation of 
the state vectors. The operator (jy^ assures us that the first path in the first column 
makes a directional change from up to right at row k. This has the effect that the 
occupation number sequence will contain at least one entry ni, where I > k. Trans- 
forming the occupation number sequence to a partition {A}, we instantly receive 
the result, {A} D {k}. 

The fact that the first path in the first column turns right at row k also means 
that the highest row that the A^th path can be when it crosses between column N 
and A^-|- 1 is fc. Thus the highest partition obtainable from lattice paths under this 
restriction are {A} = {{M)^^-^\k}. □ 



^For more details see section I.V of [8]. 



PHASE MODEL AND TODA HIERARCHY 



15 



(1 k) 

Combinatorial definitions of V'{a} • Considering the lattice path interpre- 
tation we receive, 

allowable paths in 
(M + 1) X N lattice^ 

where the lattice paths are under the condition that the first path in the first column 
makes a directional change from up to right at row k. 

Considering the lower diagonal plane partition interpretation we receive, 



wcr diag. plane part. 
N X N X M arrayt 



where the lower diagonal plane partitions are under the condition that the entry 
'Kn,! is equal to k. 

Finally, considering the ascending Young tableaux interpretation, notice that 
when we transform from the lower diagonal plane partition to the Young tableau, 
the fact that ttn,! — k means that the weight ti equals k. Since the weight ti does 
appear, as vi is not present, we can simply consider the skew partition {A — fc} to 
generate the tableaux, leading to, 

/ (l.fe) ~2t2+M -2tM+M , \M a I -2 -2\ 

(46) '^W ^2 •••«JV " =(v2...Vn) S'{A}/{fc}(«2 ,---,t^iv) 

™{ A — fc } 



3.2. Boundary correlation functions - I. Consider then the following boundary 
correlation function, 



/ T-T-iV \ M 

(47) {*iJVM)=U^F^ E 5{,}(K})5{,}/{,)(^^^...,«;^^) 

Vlij^l^j/ /Aici(ii.n(«-i).fct 



■' / {>,}c{(j\f)(«-i),fc} 



which calculates all the weighted non crossing column strict lattice paths on an 
(M + 1) X 2N lattice with the first path in the first column turning right at row k. 
Compare it to any of the wave-functions calculated earlier, and concentrate on the 
particular row, s = n — M —1 = m + N — 1, of the wave-matrix W^^"^ . If we restrict 



the variables in the usual way (eq. 40 ) , and take the vi oo limit we immediately 
obtain, 

\im r{n-M-l)w^°\n-M-l)^ ^ (K})S{,}/{fe} ({^-^j) 

{A}C{(M)(™-1) ,fc} 
{X}2{fc} 

'■prJV \M 

n,=2 J 

for < fc < M. Thus the wave-vector, given by the s = n — M— l = m + N — 1 
row of the wave-matrix, in the vi —^ oo limit gives exactly (up to a multiplicative 
factor) all the weighted non crossing column strict lattice paths on an (Af + 1) x 2N 
lattice with the first path in the first column turning right at row fc, < fc < M. 

Single determinant form for the wave-functions. When introducing the 



scalar product, we gave a single determinant form given by eq. 38 From this ex- 
pression, it is possible to obtain a single determinant form for the wave-functions 
considered above[Hl 

^•^The details below arc given in section VI of [5] to obtain single determinant expressions of 
boundary 1-point correlation functions. We expand upon these results shortly. 
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To achieve this, we first examine the operator C{v) briefly. More exphcitly, we 
are interested in the parts of C{v) that contain only (pj operators, 

M 

(48) C{v) = v'^~^^ (pj + terms that contain operators (p] 

3=0 

Thus when C{v) acts on the conjugate vacuum we obtain, 

M 

(49) {0\C{v)=v"^v-'^0\cj>^ 

We use this to express the scalar product as the following weighted hnear sum of 
correlation functions, 

M 

(50) (^M|*M)=<^^;r''(*ij^|fM) 

Therefore, if we expand the single matrix form for the scalar product as a poly- 
nomial in v\, the coefficients will reveal a single matrix form for the correlation 
functions / wave- functions. 

Using the following symmetric polynomial identity, 

(51) hp{{v'^},v^j) - hp{{v'^},vl) = - vl) hp-i{{v^},v'j,vl) 
where {v^^vf.} ^ {v^}, we can apply the following row operations, 

Rj^ - RN-k+i ,l<jk<N-k,l<k<N-l 

to completely eliminate wi from the Vandcrmondc expression in the scalar product. 
Additionally, applying the following polynomial expansion of hp{{v^},Vj), 

V 

9=0 

to all entries in the determinant which contain v\, we receive, 

M 

(52) (*m|*m) = ^r^'f^^idet 

9=0 

where. 



hM+N-l{ul, v|) 



j=2,...,N 
k = l,...,N 



(53) ^.. = u-^u-^ n n 

il = l 12=2 l<3i<fci<-;V Jl Kl 2<j2<k2<N 32 k2 

which gives a single determinant form for the (restricted) wave-functions, w^^^ 

3.3. Wave-functions - II. Considering the 'w^°°^ class of wave- functions, using 
the definitions given previously we have, 

T-{s)wi°°\s) = {-!)'' X{x}/{i''}{x)x{X}{-y) 

{A}C(n-s)(<'-'") 

where we have used the following result, 

(54) Cj{-d3)X{X}{x) = (-l)'X{A}/{in(^) 
for all partitions {A} such that {A} 3 {l-'}. 

Thus the lower triangular wave-matrix has the form. 



Constructing A'^-particle state vectors. Consider the following state vector, 

s(m)...s(u^-fc)(<Al)'|o) = |*il'>) 



PHASE MODEL AND TODA HIERARCHY 



17 



The allowable partitions of this vector are given by the following result. 
Proposition 5. 

(55) l*if')= E <Af^|A) 

{A}C{(Af)("-'=),l'=} 

Proof. Consider again the lattice path interpretation of the state vectors. The 

operator(s) (^^lA^ assure us that the last k paths, labelled jq, N — k+l<q<N, 

make directional changes from right to up at row 1, column q. Thus the largest 
occupation number sequence can be, 

{no, ni, . . . , um} = {0, fc, 0, . . . , 0, iV - fc} ^ {A} C {(M)'"^-*^', 1*} 

Also, since columns {N — k + 1, . . . , N} only contain one 0^ operator each, this 
means that only columns {1, . . . , iV — fc} can contain paths in the zeroth row. Due 
to the paths being column strict the lowest occupation number sequence is, 

{no,ni,...,nM} ^ {N - k, k,0, . . . ,0,0} ^ {A} 3 {l''} □ 

Combinatorial definitions of w^xj • Considering the lattice path interpreta- 
tion we obtain, 

r{x} — ^1 ■■■^N-k 

allowable paths in 
(M + 1} X N lattice"^ 

where the lattice paths are under the condition that the last k paths, labelled jq, 
N — k+l<q<N, make directional changes from right to up at row 1, column q, 
and only columns {1, . . . , N — k} can contain paths in the zeroth row. 

Considering the upper diagonal plane partition interpretation we obtain, 

V^{A} ^ "i • • • "iV-fe 

upper diag, plane part, 
in TV X Af X A/ arrayt 

where the upper plane partitions are under the condition that the top-right most 
k X k entries are equal to one. This obviously places restrictions on the remaining 
entries, as per the conditions of a plane partition. For example, the remaining 
{N — k) X k bottom-right entries can only either be zero or one accordingly. 

Finally, whenever we biject from the upper plane partitions to the descending 
Young tableaux, the weights tjv-fc+i = • • • = ^jv = 1, and their position in the 
tableau are exactly {In = Tj^WtN-i = i\ . . . , tAr-fc+i = T^i^}. Since these 
weights do not enter the equation, due to UN-k+i, ■ ■ ■ ,un not being present, we 
can consider the skew partition {A — l*^'} to generate the tableausr^ Thus we 
obtain, 
(56) 



M 

,(2,1'=) 2ti-M 2tM-k-M I 1 '\ c / 2 2 



Ul . . . UM-k 



UN-k) 



3.4. Boundary correlation functions - II. Consider then the boundary corre- 
lation function, 
(57) 



lj = l "'■'/ {A}C{(M)(«-fc).lfc} 



^■^Incidentally, it is at this point the reason we considered the tableaux in descending order 
becomes apparent. Had we considered ascending order we would need to invert the numbers to 
obtain the required results. 
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which calculates all the weighted non crossing column strict lattice paths on an 
{M + 1) X 2N lattice with the final k paths, labelled jq, N - k + 1 < q < N, 
turning up at row 1, column N — k + l<q<N. Additionally, only columns 
1 < q < N — k can contain paths in the zeroth row. Compare the above result with 
the s = n — M = m + N row of the wave-matrix T^(°°), restricting the variables as 
usual, and taking the UN^k+i = ■ ■ ■ = = limit, 

(58) 

lim r(n-M)4°°'(n-M) = (-!)'= ^ S^,y^^,.y{{u'})S^,y{{v-'}) 

jv-fc + i<j<jv {a}c{m("-''),i'=} 

/-r-riV-fc \ A/ 

V n,=i J 

for < fc < Af. Thus the wave- vector, given by the s — n — M — m + N row 
of the wave-matrix, in the UN-k+i — ■ ■ ■ ^ — limit gives exactly (up to 
a multiplicative factor) all the weighted non crossing column strict lattice paths 
on an (Af -I- 1) x 2N lattice with the final k paths, 1 < fc < A^, labelled jq, 
N — k + l<q<N, turning up at row 1, column N ~ k + 1 < q < N and only the 
first N — k columns can contain paths in the zeroth row. 

Single determinant form for the wave-functions. We begin by examining 
the operator B{u), as we are interested in the parts of B{u) that contain only 0| 
and operators. 



(59) 



r M M-2 

B{u) = u-'' E"''-^] + E n^^'4<^i<Pl2 

lj=0 i=0 

-I- terms that contain operators (pj , j £ {2, . . . , M} 



1-point boundary functions. Following the corresponding workings from section 
3.2, we can obtain the scalar product as the following weighted linear sum of 1-point 
boundary correlation functions, 

M 

(60) {*m|*m) = 

Explicitly expanding the determinant expression as a polynomial in we obtain, 

(61) (*m|^m^) = fifljvdet \hM+N~k{{u^}k,v']),hM-q{{u^}N-i,v'])\ j=i,...,n 

^ ' fc = l,.-.,N_i 

where {u^}k = {uf , . . . , m^} and fl^j^ is the equivalent expression of eq. 53 



2-point boundary functions. We now build upon eq. |6l] and use eq. |59] to 

consider the following quantity, 

M M-2 

(62) B(n^_i)<^llO) = uJ,'!,J2^'^-iM\0)+^nT'^ E ^n-iMI^) 

Hence we can express the 1-point correlation function as the following 

linear sum of 2-point correlation functions, 

^^'^ + E {(*M|*if ^) + (*M|*r^'°*)} 

where the coefficient of u^^_!^^ is {^mI'^m^^) + ('^m\^m'°^)- Thus in order to 
obtain {^m\^]^P), we need to &-st find (^'^l^'i/ "^). 
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This can be achieved by expanding (^'7\/|^'J^,j ) as a series in 



J=0 



Substituting q — into eq. 61 and expanding as a polynomial in u'j^_-^ we obtain, 
(64) 



k = l,...,N-2 



where Cjk = hM+N-k{{u'^}k,v'j). 

With the above result, we now expand the determinant form of (^'jv/|^'^}^), (the 
q — 1 case of eq. 61 1, as a polynomial in 

M+l 1 
91=2 92=0 

X det \cjk,hM~q2 + l{{u^}N~2,Vj),hM-qi+l{{u^}N~2,Vj)] j=i,..,,]v 

k = l,...,N~2 

where the indices (51,92) — (9,1), and (g+ 1,0), 2 < 5 < M, give us the 
sum {"ii m\^m~^'^^) + (^mI^m'"^)- Since we already have the explicit form of 
{"i' mI'^m'^^) , given in eq. 64 this leaves us with the required result. 



(65) (*iVll*M'^^) = f^fl„,ujv-ldet [Cj-fc, ftAf({li^}jV-2, l"!), ftM-9({li^}jV-2, l"!)] j = l 



k = l...-,N-2 



n-point boundary functions. Given the previous examples, we present the 
following result. 



(66) 



= 1,...,N" 
k = l,...,N — r 



Proposition 6. 

'iAf-r„_i+n-2({u^}jV-n, I'j ), . • . , flM-ri ({ } JV_„ , «| )] 3 = 

n G {0, 1, . . . , M} , ra G {0, 1} , . . . , r„ e {0, 1} 

n > 7-2 > ■ ■ ■ > r„ , 1 < n < N 

Com men t. The above statement can be proven using induction. By assuming 
that eq 



66 



is true, we obtain (^'mI^'m ' ),0<q<n — 1, as the following 



weighted sum of (n + l)-point correlation functions 

A/-1 



(*A/i* 



A/|*M 



{M,1^-1,01). -M 



(67) 



J = l 

Additionally for {^M\'^fP), {q = n), we have, 

M 

(68) (*A/i*ir>) = E^^^'-n^*"!*!/""') 



j=0 



In order to verify the proposed result we need to derive, using polynomial expansion 
method(s) on eq. 66 the explicit determinant forms for the following expressions 
(referred to as step 1, 2 and 3), 

• Step 1) {"^mI^m'^ ^), the coefficient oi u^^^^-' in eq. 



. Step 2) (vI/M|v|/ir"'^°'^'>) 



and (^'jul* 



and 



'^N-n in eq. 



67 



Ml'i'M 



{M,i"^'',a''}\ 



68 



, the coefficients of u 



-M 
N-n 
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• Step 3) (^-mI^'^" '''°''^) and (^'mI^'m^^'^" ' ''°'^'^), the coefficients of 



-M+2j 

in eq. 



67 



It is necessary in the third step to determine which of the determinant expressions 
corresponds to which correlation function. This is achieved through considering 
values of q where one term is already known from a previous result. As an example, 

consider q = n — 1, where we receive the sum (^'a/I^'^^'^'" ^) + ^). 

In this case we have already obtained (^'mI^'^/^^'" ^) from the first step. For 

q = n — 2 we receive the sum (\1/m|^'m^ ^) + m\^m^^'^'° where we 

know the value of (^'m l^i/^^'^'" '^), 1 < j < M — 2, from the previous (g = n — 1) 
calculation, and we know j — M — 1 from the second step. Carefully following this 
argument for all values of q we complete the proof by induction. 

Thus using the results of eq. |66[ we obtain the single determinant form for the 
restricted wave-functions, w^°°\ 

4. Discussion 

The main result of this work is the correspondence between the 2-Toda wave- 
functions and the boundary correlation functions of the phase model. The weighted 
sum of the wave-functions analyzed in this work can be thought of as the action of 
a single vertex operator [10] on the 2-Toda r-function, 

Txi\)r{s,x,y) _ ^fc -.(oo)^„^ ^yW^js + l,x,y) _ -.(O) 



r{s,x,y) ' r{s,x,y) 

where T^/y{X) = exp|- }. It is a pertinent question as to 

whether there is a combinatorial interpretation of a r-function that has been acted 
on by more than one vertex operator. Alternatively, it is perfectly natural to 
speculate whether non boundary correlation functions of the phase model have 
natural correspondences with fundamental objects of the 2-Toda hierarchy. It is 
the author's intentions to examine these and related questions in the future. 

Acknowledgements 

The author would like to thank O Foda for discussions, and M Wheeler for 
making him aware of the plane partition/tableau bijection, without which this 
work would not exist. This work was supported by the Dept. of Mathematics and 
Statistics, The University of Melbourne. 

References 

[1] O Foda, M Wheeler and M Zuparic, J. Stat. Mech. (2009) P03017 

[2] O Foda, M Wheeler and M Zuparic, Nuclear Phys. B 820 [FS] (2009) 649-663 

[3] V E Korepin, N M Bogoliubov and A G Izergin, Quantum inverse scattering method and 

correlation functions Cambridge University Press (1993) 
[4] N Bogoliubov, A Izergin and N Kitanine, Nuclear Phys. B 516 (1998) 501-528 
[5] N Bogoliubov, J. Phys. A: Math. Gen. 38 (2005) 9415-9430 

[6] K Ueno and K Takasaki, Advanced Studies in Pure Mathematics 4 (1984) 1-95 
[7] K Takasaki, Advanced Studies in Pure Mathematics 4 (1984) 139-163 

[8] I G Macdonald, Symmetric functions and Hall polynomials Oxford University Press (1995) 
[9] N V Tsilevich, Fund. Anal. Appl. 40 No. 3 (2006) 207-217 
[10] T Miwa, M Jimbo and E Date, Solitons Cambridge University Press (2000) 

Department of Mathematics and Statistics, University of Melbourne, Parkville, Vic- 
toria 3010, Australia. 

E-mail address: mzup@ms.uiiimelb.edu.au 



